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ABSTRACT 

We propose a new, simple but powerful algorithm to analyze the gamma-ray 
burst temporal structures based on identifying non-statistical variations 
("peaks") in the time histories. Detailed analyses of the bursts from the third 
BATSE catalog show that ~ 30 bursts have more than 20 peaks individually. 
Upon identifying most of the peaks in those bursts, we show that the peak 
fluence Si and peak interval Si distributions within each burst are consistent 
with log-normal distributions. Furthermore, we show that Gaussian (in linear 
space) and power-law distributions for peak fluences are ruled out, as is the 
Poisson distribution for peak intervals. 

Subject headings: gamma rays: bursts 
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1. Introduction 

Though the origin of Gamma-ray Bursts (GRBs) remains elusive, the abundant 
data collected by BATSE on Compton Gamma-ray Observatory (CGRO) have provided 
many constraints on the physical modeling of these events. Thanks to BATSE's high 
temporal resolution and high signal-to-noise ratio, temporal analyses of bursts have began 



to bear fruit. Notably, the bi- modal distribution in the burst durations (jKlebesadel 



1992| ; |Kouveliotou et al. 1993| ); the claimed time-dilation of peak width between strong and 



weak bursts ( [IN orris et al. 199~3| although see |Mitrofanov et al. 1994j ) ; and that the peak 



width narrows with energy following a power-law ( Fenimore et al. 1995| ). Though the exact 



physics is yet unknown to satisfy all these observations, they are important probes into the 
burst emission processes. 

Here in this Letter, we propose a new, simple but powerful peak-finding-algorithm 
(PEA) that can identify the "peaks" (or non-statistical variations) within a GRB. Using 
this PFA, we uncover that the distributions of peak fluences and intervals within each burst 
are log-normal and several tests are performed to confirm our findings. 



2. The Peak-Finding Algorithm 

GRB time histories show vast diversities and Figure [I] shows portions of two complex 
bursts: trigger 678 and 1606. Norris et al. (1996) have used a pulse-fitting algorithm based 
on x 2 -fitting of individual pulses. However, their method fails to converge when the burst 
is very complex, particularly for the two bursts shown here. 

The detailed procedure to find all the peaks and valleys in a burst is as follows: (a) 
First, we fit the burst background using a linear function B{t) to the pre- and post-burst 
regions, (b) During the burst, every count bin that has more counts than the neighboring 
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bins (both sides) is a candidate peak with count C p at time t p . (c) We then search on both 
sides of each candidate peak for counts C\ (at t{) and C 2 (at t 2 ) so that the conditions 
C p — Ci 5 2 > N vaI yJC p ~ are satisfied, (d) The search will stop either when both C\ and C 2 
are found, in which case C p becomes a "true" peak, or when counts higher than C p (on 
either side of t p ) is encountered, in which case C p is not a true peak and discarded. After 
this step, all the peaks (iV fc in total) should have been identified, (e) We then locate the 
minima between two successive peaks as valleys. Note that C% and C 2 are not necessarily 
the valleys. The valley at the beginning (end) of a burst is chosen from the location where 
counts start deviating from (dimming into) the background. 

The above procedure has been implemented and in Figure III we present the selected 
peaks (up-triangles) and valleys (down-triangles) for burst 678 and 1606. Note that only a 
small portion of the time histories is shown for each burst and there 58 and 35 peaks in 
total (with iV var = 5) for burst 678 and 1606 respectively. 

Assume that the peak is C p (t p ) and two neighboring valleys are Cx(tx) and C 2 (t 2 ) after 
applying the peak-finding algorithm, as schematically shown in Figure |2| along with the 
background level B(t). The interval between adjacent peaks, or the waiting time between 
successive peaks, is 5i = t Pi+1 — t Pi with i = 1, . . . , iV fe — 1; and the count fluence within 
the zth peak is defined as Si = J2tl [C(t) ~ B{t)]. Note that we use count fluence rather 
than photon energy fluence, though they are simply related by the mean photon energy if 
it remains approximately the same throughout the burst. So, there are N^ — l intervals but 
Nk peak fluences in each burst. 

The number and position of peaks within a GRB, as one might suspect, depends 
somewhat on the exact value of iV var , but the effect is small when 3 < iV var < 5. In this 
Letter, we choose iV var to be 5. Tests using various iV var show that our conclusions are not 
affected by this choice. 
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3. Log-Normal Distributions of Peak Fluences and Intervals 

A two-parameter log- normal distribution is represented as follows ( Aitchison fc Brown 



19571; Crow & Shimizu 19881): 




X < 

where f(x) is the probability density function for logx, fi and a 2 are the two parameters 
corresponding to the sample mean and variance of {logXj}. The cumulative log- 
normal distribution is | erfc(y), where erfc(y) is the complementary error function and 
y = (n- \ogx i )/V2a 2 ~. 

Utilizing the BATSE DISCSC data (64 ms temporal resolution and 4-channel spectral 
resolution), we analyze most of the bursts in the third BATSE catalog ( [Meegan et al.| 
|1996|) and identify their peaks using our PFA with N var = 5. Since we require that each 
burst must have more than 20 peaks so that a distribution can be built, and all the peaks 
must be 5 sigma above background, 32 bursts are selected that meet the criteria. The 
requirement of more than 20 peaks in each burst is somewhat arbitrary and it is mainly due 
to the concern that the statistical significance on fewer than 20 points in a distribution is, 
in general, questionable. We have done tests by relaxing this requirement and found that 
our conclusions remain intact. 

All 32 selected bursts tend to be bright, long, complex and, on average, each of 
them has ~ 35 peaks. In the following, we concentrate on the peak fluences Si and 
peak intervals Si distributions. For a burst with A" peaks, we define two sample means 
for fluences: fii n = J2± ^og(Si)/N and jj n = J2i Si/N; and two sample variances: 
v? n = Ef (log(^) - fii n ) 2 /(N - 1) and o% = Ei (Si - Vn) 2 /(N - 1). Here, the subscript 
In represents log and n for linear. Sample means and variances for {Si} can be defined 
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similarly. 

We find that, for all 32 bursts, the distributions of the peak fluences {Si} and the 
peak intervals {5i} in individual bursts are consistent with log-normal distributions (i.e. the 
differential distributions of {log(S'j)} an d {log(<5j)} are consistent with Gaussian). We reach 
this conclusion by performing x 2 -fitting. The peak fluences and intervals from each bursts 
are binned (5 bins) and fitted with log-normal distributions with two parameters (mean and 
variance). In Figure |3|a and 6, we give the x 2 -fitting probabilities for {Si} and {<5j} for each 
burst, respectively. Note that the probabilities are reasonably uniformly distributed between 
and 1, implying that log-normal is, in general, an acceptable hypothesis. However, we 
caution that x 2_ fitting might not be an accurate measure of statistics here due to the small 
number of points in each bin though they do satisfy the limit given by Lindgren (1976) that 
the sample size should be roughly 4 or 5 times the number of bins. 

Perhaps the peak fluence distribution is similar to log-normal, but instead follows some 
other distributions, such as a truncated power-law. We have performed fitting to peak 
fluences assuming they are distributed as power-laws (again with two parameters, the slope 
of the power-law and the lower cutoff). Their probabilities are shown in Figure (notice 
that probability is now in logarithmic scale). More than half of the bursts have probabilities 
less than 10 -3 and the acceptable probabilities for some bursts may just be due to their 
small number of peaks (note the trend for probability becoming exceedingly smaller when 
peak number increases). 

Even though peak fluence and peak interval distributions in individual burst are 
consistent with log-normal distribution, their confidence is nevertheless affected by the 
small number of peaks. A more stringent test on peak fluence distribution can be 
performed as follows: for each burst, we rescale each Si so that the new fii n = \i n = 
and af n = a\ = 1. As a result, all 32 bursts now have the same ji{= 0) and cr 2 (= 1) in 
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the peak fluence distribution. We then put all the peaks from all 32 bursts together and 
sort them in the ascending order for the scaled Si. This combination of all peaks from 
all bursts after the rescaling (1107 peaks total) removes the uncertainty in small number 
of peaks. In Figure and b, we present the cumulative distribution of rescaled log Si 
and Si as compared to the cumulative Gaussian distribution with mean = and variance 
= (iV pk - N hur )/(N pk - 1) « 0.97, where iV pk = 1107,iV bur = 32 are the total number of 
peaks and bursts, respectively. This extremely good fit in Figure |]a strongly indicates that 
the peak fluences in GRBs are distributed log-normally rather than normally (Figure |4]&) . 

The same procedure can be applied to the peak interval distribution and the results 
for log 5, are shown in Figure [|c, it is well fitted by the log-normal distribution. Special 
attention should be taken, however, in making the comparison with the hypothesis that 
peak intervals are distributed randomly. The probability for having another peak after time 
Si is e -^l<^> (i.e. the Poisson probability for zero event), and the corresponding cumulative 
distribution is 1 — e~ x , where < x < oo. Therefore, we rescale the peak intervals in each 
burst by its mean < Si > and put all the peaks from all bursts together and plot it against 
1 — e~ x . This is shown in Figure |4]<i. The inconsistency is obvious, therefore, the hypothesis 
that peaks appear randomly can be ruled out. 

Though the fitting to log-normal distribution is extremely good, small deviations in 
Figure |4]a and |]c are still visible. Possible explanations include the following: For peak 
fluence, our definition of Si (the area enclosed by two successive valleys) may be a slight 
underestimation, especially for large peaks due to the fluences lost in the tails which extend 
into the neighboring peaks. However, this effect is mutual for adjacent peaks so the net 
difference can be estimated to be very small (a few percent at most). This is enough 
to account for the tiny asymmetry (bright end) in peak fluence distribution. The small 
asymmetry in peak interval distribution, we believe, is caused by the limited temporal 
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resolution (64ms) of the data. In fact, the peak interval has to be larger than a few bins 
(by definition) whereas equation [I] assumes x can be arbitrarily close to 0. More precisely, 
a three-parameter log-normal distribution ([Aitchison fc Brown 1957| ) can be employed but 
we omit any detailed comparison here as the real correction can only come from analyzing 
the high temporal resolution data (such as Time- Tagged Events (TTE) data of BATSE). 

McBreen et al. (1994) have also suggested that the intervals between peaks in GRBs 
may be distributed log-normally, based on their analyses of 34 long, intense KONUS bursts 
with 0.25 s time resolution (1 s resolution after 34 s in a burst). Their analyses and 
results, however, are very different from ours in several aspects. By using a rigorous peak 
selection process, we have identified individual entities within the GRB time history that 
are distinct, and found that each individual burst has a log-normal distribution. In contrast 
to our ~ 35 peaks per bursts, McBreen et al. (1994) had only 4-5 "peaks" on average. 
More questionably, peak intervals from different bursts with different brightnesses were put 
together to form one distribution with no normalization. We find that great caution has 
to be taken when putting together quantities from different bursts due to the bias that 
fewer peaks on average will be selected from weaker bursts. No alternative, non-lognormal 
distributions were considered in their analyses either. 



4. Summary 

Our findings that peak fluences and peak intervals in individual GRB are distributed 
log-normally probably have far-reaching implications on the production of GRBs and 
consequently their origin, though it is beyond the scope of this Letter. The non-random 
appearance of peaks seems directly pointing to a central engine, and interestingly, similar 
behavior has also been found in the interval (waiting time) distribution for the successive 
events from the soft gamma-ray repeaters (SGRs) ( [Hurley et al. 1994]) , which have been 
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suggested to come from the neutron star quakes flCheng et al. 1996| ). Further analyses on 
the statistical properties of GRB time histories, such as the detailed comparisons with those 
found in SGRs, will be presented in the future publications. 

Useful discussions with B. Cheng, S. Colgate and R. Epstein are appreciated. H.L. 
gratefully acknowledges the support by the Director's Postdoc Fellowship at Los Alamos 
National Lab. This work is partially supported by the GRO Guest investigator program. 
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Fig. 1. — Portions of time histories for 2 complex GRBs 678 and 1606. Filled up-triangles 
show chosen peaks and down-triangles show valleys using our peak-finding algorithm. There 
are 58 and 35 peaks in total for burst 678 and 1606, respectively, although only part of the 
time histories are shown. 

Fig. 2. — A schematic presentation of a selected peak C p (t p ) and two neighboring valleys 
Ci(ii) and 6*2(^2) in a burst time history. The condition C p —C\^ > A" var ^Jc~ p with N vax > 3 is 
satisfied, ensuring C p is a true peak not a random statistical fluctuation. The line -Bi, p ,2 shows 
the background level. The peak fluence is defined approximately as Si = Y?t\ — B{t)} 
and peak interval 5, = t Pi+1 — t Pi . 

Fig. 3. — The x 2 probability (filled dots) for each burst's peak fluences (a) and peak intervals 
(b) to be log- normally distributed versus the number of peaks in that burst. Also shown is 
the probability (triangles) if peak fluences in each burst are distributed as power-laws (c). 
Note the logarithmic scale of probability in plot c. Each dot represents one burst (32 total). 
The consistency with log-normal for fluence and interval distributions and the inconsistency 
with power-law for peak fluences are apparent. 

Fig. 4. — The rescaled and sorted peak fluences and peak intervals from all 32 bursts (dots in 
all panels) are compared to both the log-normal distribution (dashed line in a, hardly visible, 
and c) , the Gaussian distribution (dashed line in b) and Poisson interval distribution (dashed 
line in d). The rescaling is done by (x— < Xi >)/ < (x— < X{ >) 2 > x / 2 , whereas Xi is 
log(Si) in plot a, Si in b, and log(5j) in c, respectively. It is clear that burst distributions are 
consistent with log-normal whereas Gaussian hypothesis for fluences and Poisson hypothesis 
for intervals are ruled out. 
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